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We point out that infiaton decays can be a copious source of stable or long-lived particles \ with 
mass exceeding the reheat temperature Tr but less than half the infiaton mass. Once higher order 
processes are included, this statement is true for any \ particle with renormalizable (gauge or 
Yukawa) interactions. This contribution to the \ density often exceeds the contribution from thermal 
X production, leading to significantly stronger constraints on model parameters than those resulting 
from thermal \ production alone, particularly in models containing stable charged particles. 
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According to inflationary models (lj, which were 
first considered to address the flatness, isotropy, and 
monopole problems of the hot Big Bang model, the Uni- 
verse has evolved through several stages. During infla- 
tion, the energy density of the Universe is dominated 
by the potential energy of the infiaton and the Universe 
experiences a period of superluminal expansion. Immedi- 
ately after inflation, coherent oscillations of the infiaton 
dominate the energy density of the Universe. These os- 
cillations eventually decay, and their energy density is 
transferred to relativistic particles; this reheating stage 
results in a radiation-dominated Fricdmann-Robcrtson- 
Walker (FRW) Universe, as in the hot Big Bang model. 

Initially reheating was treated as the perturbative, one 
particle decay of the infiaton with decay rate Td, result- 
ing in Tr ~ (r^A/p) 1 / 2 for the reheat temperature |p^|, 
where Mp = 2.4 x 10 18 GcV is the reduced Planck mass. 
Tr should be low enough so that the original monopole 
problem is avoided. Moreover, in many supcrsymmetric 
models Tr < 10 7 — 10 9 GcV, in order to avoid gravitino 
overproduction which would destroy the success of nucle- 
osynthesis B. Later it has been noticed that the initial 
stages of infiaton decay might involve non-perturbativc 
resonance processes Q. They typically lead to a highly 
non-thermal distribution of particles, including inflatons 
with large momentum 0]. However, after sufficient red- 
shifting the energy density of the Universe would again 
be dominated by non-relativistic, massive particles. It 
is therefore generally believed that an epoch of (pertur- 
bative) reheating from the decay of massive particles (or 
coherent field oscillations, which amounts to the same 
thing) is an essential ingredient of any potentially realis- 
tic cosmological model Q . In what follows we generically 
call the decaying particle the "infiaton" , since we are (al- 
most) sure that inflatons indeed exist. Note also that in a 
large class of well-motivated models, where the infiaton 
resides in a "hidden sector" of a supergravity theory , 
its couplings are suppressed by inverse powers of Mp, and 
hence are so weak that infiaton decays are purely pertur- 
bative. However, it should be clear that our results hold 
equally well for any other (late) decaying particle. 



Even before all inflatons decay, their decay prod- 
ucts form a plasma which, upon a very quick ther- 
malization, has the instantaneous temperature H T ~ 
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^g* ' "HT^MpJ , where H is the Hubble parameter 

and denotes the number of relativistic degrees of free- 
dom in the plasma. This temperature reaches its maxi- 
mum T max soon after the infiaton field <p starts to oscil- 
late, which happens for a Hubble parameter Hi < m^, 
with ?7!.0 being the frequency of infiaton oscillations about 
the global minimum of the potential. We will assume that 
all infiaton decays can be described by perturbation the- 
ory in a trivial vacuum, which implies T max < to^ /2. (The 
resulting upper bound on also implies that a vacuum 
expectation value of the infiaton field does not induce 
large masses to the particles to which it couples.) How- 
ever, T ma x can be much larger than Tr. As long as 
T > Tr the energy density of the Universe is still domi- 
nated by the (non-relativistic) inflatons that haven't de- 
cayed yet. The Universe remains in this phase as long 
as H > Td- During that epoch particles x with mass 
Tmax > fn x > can produced copiously from the 
thermal plasma [pldTT). Here we point out that \ par- 
ticles can also be produced directly in infiaton decays. 
We will show that the x abundance from infiaton decay 
often exceeds that from thermal production, even if the 
branching ratio for <fi — ► \ decays is very small. 

We begin our argument by pointing out that T max is 
frequently well below m^. This is important, since ther- 
mal production is obviously only efficient if m x < T max , 
while infiaton decay can produce pairs of x particles as 
long as m x < m^/2. For perturbative infiaton decay 
thcrmalization increases the number density and reduces 
the mean energy of the decay products. Complete thcr- 
malization (i.e. both chemical and kinetic) therefore re- 
quires 2 — > N reactions, which change the number of 
particles, to be in equilibrium. Since the rate for higher 
order processes is suppressed by powers of the relevant 
coupling constant a, the most important reactions are 
those with N = 3. These reactions have recently been 
studied in Ref. (L2| where the scattering of two matter 
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fermions with energy ~ (from inflaton decay) to 

two fermions, plus one gauge boson with typical energy 
E <C Tn^ , is considered. The rate for these reactions 
can be large due to the t— channel pole of the scattering 
matrix element, regulated by a cut-off on the exchanged 
momentum, naturally taken to be the inverse of the aver- 
age separation between two particles in the plasma |[L2| . 
It turns out that the largest possible T max is given by |13|| 
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Even if is near its upper bound of ~ 10 13 GeV Q, 
for a chaotic inflation model, and Tr is around 10 9 GeV 
(saturating the gravitino bound) T max will exceed Tr if 
the coupling a 3 > 10 -8 . This is easily accommodated 
for particles with gauge interactions. On the other hand, 
recall that T max < Together with eq.(|l|), taking 

a < 0.1, this gives T max < 10 11 (10 5 ) GeV for T R = 
10 9 (1) GeV. This implies in particular that there will be 
no "wimpzilla" production |S|] from thermalized inflaton 
decay products, since in this case m x > T max . 

On the other hand, for m x < 20Tr the standard calcu- 
lation |Q of the density of stable relics applies. Scenarios 
with T max > m x > 20Tr have only been investigated rela- 
tively recently in refs. ]^-[Tl|l, which studied x production 
from the thermal plasma with T > Tr. If the x density 
was always well below the equilibrium density, one finds 
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Here Vl x is the x mass density in units of the criti- 
cal density and h is the Hubble constant in units of 
100 km/(s-Mpc). We have taken the cross section for 
X pair production or annihilation to be a ~ ar./m^.. 
Note that tt x is only suppressed by (Tr/to x ) 7 rather 
than by exp(— m x /T^). A stable particle with mass 
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m x <~*j 2000 Tr • a x might thus act as the Dark Matter 
in the Universe (i.e. fl x ~ 0.3). However, eq.(Q) with 
a = 0.05 implies that T max > IOOOTr is only possible if 
T R < 2 • 10- 12 A/ P . Eq.@ is no longer applicable [§ if 
the coupling a x is so large that x reached chemical equi- 
librium; however, it can then still be used as an upper 
bound on r^ herm . 

We now discuss the direct production of x particles in 
inflaton decay. (Other mechanisms for nonthermal pro- 
duction of superheavy particles have been discussed in 
Jl5|.) Most inflatons decay at T ~ Tr; moreover, the 
density of x particles produced in earlier inflaton de- 
cays will be greatly diluted. Since inflaton decay con- 
serves energy, the density of inflatons can be estimated 
as ~ O.3g»T^/m0. Let us denote the average num- 
ber of x particles which are produced in each <fi decay by 
B{<t> — ► x)- We translate the x density at T = Tr into 
the present x relic density using the relation H 
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The x density from <j) decay is therefore fluj : 
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Eq. (Q) holds if the x annihilation rate is smaller than the 
Hubble expansion rate at T ~ Tr, which requires 
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This condition will be satisfied in chaotic inflation mod- 
els with rrifj, ~ 10 _5 Afp, if m x is large enough to avoid 
overclosure from thermal x production alone. It might 
be violated in models with light inflaton. In that case 
the true x density at Tr can be estimated by equating 
the annihilation rate with the expansion rate: 
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This maximal density violates the overclosure constraint 
il x < 1 badly for the kind of weakly interacting (a x <0.1), 
massive (m x 3> Tr and m x > 1 TeV) particles we are 
interested in. [Eq.(^) describes the maximal x density if 
X decouples at T ~ Tr. It is not applicable to WIMPs 
decoupling at T < Tr.] For the remainder of this article 
we will therefore estimate the x density from inflaton 
decay using eq.(^). 

Our remaining task is to estimate B((f> — > x)- This 
quantity is obviously model dependent, so we have to 
investigate several scenarios. The first, important spe- 
cial case is where x IS the lightest supersymmetric par- 
ticle (LSP). If m^y is large compared to typical visible- 
sector superparticle masses, <fi will decay into particles 
and superparticles with approximately equal probabil- 
ity. (This statement is true so long as the superpoten- 
tial is quadratic or higher in the inflaton superfield ]l3|.) 
Moreover, all superparticles will decay into one x parti- 
cle and some standard particle (s) at a time scale which 
is shorter than the superparticle annihilation time scale 
p6| , as long as m x > Tr, even if a x ~ 0.1. As a result, if 
X is the LSP, then B(<f> — > x) — 1j independently of the 
nature of the LSP. 

Another possibility is that the inflaton couples to all 
particles with more or less equal strength, e.g. through 
non-renormalizable interactions. In that case one ex- 
pects B{<f> — > x) ~ ~ 1/200. However, even if <j> has 
no direct couplings to Xi the rate (Q) can be large. The 
key observation is that x can be produced in (f> decays 
that occur in higher order in perturbation theory when- 
ever x can be produced from annihilation of particles in 
the thermal plasma. In most realistic cases, (j) — > f fxx 
decays will be possible if x has electroweak gauge in- 
teractions, where / stands for some gauge non-singlet 
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with tree-level coupling to cj). A diagram contributing 
to this decay is shown in Fig. 1. Note that the part of 
the diagram describing XX production is identical to the 
diagram describing XX f f transitions. This leads to 
the following estimate: 
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where C4 is a multiplicity (color) factor. The phase space 
factors have been written in a fashion that reproduces the 
correct behavior for m x — > m^/2 as well as for m x — > 0. 
Occasionally one has to go to even higher order in per- 
turbation theory to produce x particles from <j) decays. 
For example, if x has on ly strong interactions but <fi only 
couples to SU (3) singlets, XX pairs can only be produced 
in six body final states, <j) — > ffqqxx- A representative 
diagram can be obtained from the one shown in Fig. 1 
by replacing the x lines by quark lines, attaching an ad- 
ditional virtual gluon to one of the quarks which finally 
splits into XX- The branching ratio for such six body 
decays can be estimated as 
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Another example where XX pairs can only be produced 
in <j> decays into six body final states occurs if the inflaton 
only couples to fields that are singlets under the standard 
model gauge group, e.g. right-handed (s)neutrinos vr 
jL7| . [Since vr decays very quickly, B(<j> — > vr) ~ 1 does 
not cause any problem.] Since vr only has Yukawa inter- 
actions, the factor a 2 ^ in eq.(||) would have to be replaced 
by the combination of Yukawa couplings A 2 H A 2 /(167r 2 ). 
If 2m x < m VR , XX pairs can already be produced in four 
body final states from Vr decay. The effective <f> — > x 
branching ratio would then again be given by eq. (0) , with 
m^ replaced by m VR in the kinematical factors. 

Finally, in supergravity models there in general ex- 
ists a coupling between <f> and either x itself or, for 
fcrmionic to its scalar superpartner, of the form 
a (m ( f > m x /Mp) (f>xx + h.c. in the scalar potential jfjj. A 




Fig. 



1: Sample diagram for x production in four-body 
inflaton decay. 



reasonable estimate for the coupling strength is |TJ] a ~ 
(4>)/Mp, unless an R— symmetry suppresses a. Assuming 



that most inflatons decay into other channels, so that 
Td ~ \fShT-g_jM-p remains valid, this gives 



B(cj> x) 
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The production of x particles from inflaton decay will 
be important for large m x and large ratio ?ti x /Tr, but 
tends to become less relevant for large ratio m^/m x . 
Even if m x < T max , x production from the thermal 
plasma (0) will be subdominant if 
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The first factor on the r.h.s. of (|l(J) must be < 10~ 6 in 
order to avoid over-production of x from thermal sources 
alone. Even if <fi —> x decays only occur in higher orders 
of perturbation theory, the l.h.s. of (1^) will be of order 
10~ 4 (10~ 10 ) for four (six) body final states, see eqs.(0), 
(g); if 4> — > XX decays at tree-level, the l.h.s. of ( [To| ) will 
usually be bigger than unity. We thus see that even for 
10 13 GeV, as in chaotic inflation models, and for 
10 3 Tr, x production from decay will dominate if 
m x >10 7 (10 10 ) GeV for four (six) body final states. As a 
second example, consider LSP production in models with 
very low reheat temperature. The LSP mass should lie 
within a factor of five or so of 200 GeV. Recall that in 
this case B{<j) — > x) — !• Taking a x ~ 0.01, we see that 
X production from decay will dominate over production 
from the thermal plasma if m^ < 6 • 10 7 GeV for Tr = 1 
GeV; this statement will be true for all m<f, < 10 13 GeV 
if T R < 100 MeV. 

Let us now assume that eq.(^) indeed gives the domi- 
nant contribution to x production in the early Universe, 
and investigate the resulting constraints on model pa- 
rameters. As well known, any stable particle must sat- 
isfy Vl x h 2 < 1, since otherwise it would "overdose" the 
Universe. For example, in case of a neutral LSP with 
m x ~ 200 GeV, eq.® with B(<j> -> x) = 1 implies 
fn<f>/Tp_ > A ■ 10 10 . Such a large ratio to^/Tr in turn re- 



quires Td < 10 



l/Mp, which indicates that <p would 



have to decay through higher dimensional operators. Of 
course, this constraint is no longer valid if x reaches equi- 
librium with the plasma at temperatures < Tr . 

Another Dark Matter candidate is a very massive par- 
ticle, with m x ~ 10 12 GeV; decays of this particle could 
give rise to the observed very energetic cosmic rays ]l9| if 
their lifetime is > 10 8 times the age of the Universe. We 
noted above that such massive particles cannot be pro- 
duced thermally in any realistic model of inflation. On 
the other hand, eq.(^) shows that inflaton decays might 
very easily produce too many of such particles. Taking 
m$ = I0m x — 10 13 GeV, we see that we need a branch- 
ing ratio as small as 5-10 -8 GeV/TR, which implies quite 
a severe upper bound on Tr even if x pairs can only be 
produced in six body decays of the inflaton. Even taking 
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Tr = 1 MeV, the lowest value compatible with successful 
nucleosynthesis, this requires B((f> — > x) < 10~ 4 . Finally, 
if x is produced only through Mp suppressed interac- 
tions, eq.(|) implies a 2 < 3.5 • 10" 6 GeV • M v T K /m%, 



which again gives a very tight constraint if m x 
GeV. 

In some cases other considerations 
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give an even 
stronger constraint on Cl x . For example, the abun- 
dance of charged stable particles is severely constrained 
from searches for exotic isotopes in sea water pj| , e.g. 
n x h 2 < 1(T 20 for 100 GeV<m x <10 TeV; for heavier par- 
ticles this bound becomes weaker. This bound imposes 
very severe constraints on supersymmetric models with 
stable charged LSP. Fixing again m x = 200 GeV from 
considerations of naturalness, to^/Tr > 4-10 30 B(cf) — > x) 
is required. This is clearly incompatible with the limits 
Tr > 1 MeV, 77i0 < 10 13 GeV, even if <\> — > x decays re- 
quire six body final states, see eq.(||). We saw above 
that arranging x to have been in equilibrium at Tr does 
not help. Finally, the relic density of charged LSPs that 
were in thermal equilibrium at T < Tr is too large by 
more than ten orders of magnitude. Eq.([|) shows that 
the situation for larger m x would be even worse. We 
thus conclude that in models where at least a significant 
fraction of the present entropy of the Universe originates 
from inflaton decay, a stable charged LSP can only lead 
to an acceptable cosmology if it is too massive to be pro- 
duced in inflaton decays. 

Our calculation is also applicable to entropy-producing 
particle decays that might occur at very late times. If x is 
lighter than this additional <j>' particle |[T6| , all our expres- 
sions go through with the obvious replacement <f> — > <f>' 
everywhere. More generally our result holds if <f> de- 
cays result in a radiation dominated era with Tr > my . 
If <j>' is sufficiently long-lived, the Universe will eventu- 
ally enter a second matter-dominated epoch. <j>' decays 
then give rise to a second epoch of reheating, leading 
to a radiation-dominated Universe with final reheating 
temperature Ta f , and increasing the entropy by a factor 
to^'/Trj. This could be incorporated into eq.(|J) by re- 
placing Tr — > TpJTp^/m^' > Trj. Our result regarding 
a stable charged LSP would remain valid in such a sce- 
nario even if m x > , since the lower bound of ~ 1 
MeV which we used now applies to Tr ( . . The only way 
out would be to allow <j)' to be essentially the only decay 
product of (j), where <j) itself does not have renormalizable 
interactions with standard particles and their superpart- 
ners (so that higher order (f> decays are negligible) and 
2m x > m^i . However, there is presently no motivation 
for considering such baroque models. 
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